Anderson localization of phonons in dimension d = 1,2,3 : 
finite-size properties of the Inverse Participation Ratios of eigenstates 
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We study by exact diagonalization the localization properties of phonons in mass-disordered har- 
monic crystals of dimension d = 1, 2, 3. We focus on the behavior of the typical Inverse Participa- 
tion Ratio ¥2(1*], L) as a function of the frequency uj and of the linear length L of the disordered 
samples. In dimensions d = 1 and d = 2, we find that the low-frequency part uj — > of the spec- 
trum satisfies the following finite-size scaling LYi(uj,V) = Fd=i{L 1 ' 2 u>) in dimension d — 1 and 
L 2 Y2(oj, L) = i 7 d = 2((lnL) 1/ ' 2 a)) in dimension d = 2, with the following conclusions (i) an eigenstate 
of any fixed frequency uj becomes localized in the limit L — > +00 (ii) a given disordered sample of 
size L contains a number Ndeioc(L) of delocalized states growing as N deloc (L) ~ L 1/2 in d = 1 
and as Ndeioc(L) ~ L 2 /(lnL) in d — 2. In dimension d — 3, we find a localization-delocalization 
transition at some finite critical frequency uj c (W) > (that depends on the disorder strength W). 
Our data are compatible with the finite-size scaling L d{2) Y 2 (uj,L) = F d=z {L 1/u {uj - w c )) with the 
values D(2) ~ 1.3 and v ~ 1.57 corresponding to the universality class of the localization transition 
£NJ ' for the Anderson tight-binding electronic model in dimension d = 3. 
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INTRODUCTION 



Since its discovery fifty years ago [l], Q Anderson localization has remained a very active field of research (se e 
for instance the reviews and has been recently realized in experiments with atomic matter waves [Tol . flll j . 

According to the scaling theory pjj], there is no delocalized phase in dimensions d = 1,2, whereas there exists a 
localization/delocalization at finite disorder in dimension d > 2. However the notion of Anderson localization is 
not limited to quantum electrons models, but also applies to classical waves in disordered media [3, [l3| including 
*^ . acoustic waves (see for instance an d references therein), electromagnetic waves (see for instance [l8l [l9j and 

references therein), and hydrodynamical waves [20]. Among the classical disordered models that are expected to 
display Anderson localization, the oldest problem is the phonon problem in the presence of random masses coupled 
by fixed spring constants, which has been introduced by Dyson [2l| even before Anderson's paper [TJ. After studies 
concerning the one-dimensional case (see the review [22j and references therein), an analysis of disordered elastic 
media via a non-linear sigma-model [23| has predicted results similar to the scaling theory of Anderson localization 
[l2j : all finite-frequency phonons are localized in dimension d < 2, whereas there exists a finite critical frequency 
uj c > in dimension d > 2 that separates delocalized modes uj < uj c from localized modes uj > u> c . However, in 
contrast to electron models where many numerical studies have checked in detail these predictions and more refined 
properties like multifractal properties at criticality (see the review H), the same effort to characterize the statistics of 
eigenstates has not been done for the phonon problem. In particular in dimension d — 3, the numerical studies we are 
aware of find that almost all states are delocalized, whereas localized states appear only near band-edges [23 - |27| . In 
addition, the universality class of the transition does not seem completely clear : the reported numerical values of the 
critical exponents are sometimes the same as for the Anderson electron transition [24| . but are sometimes different 
[2ol [26} . The aim of the present paper is thus to revisit the problem of phonon localization in dimension d = 1,2,3 
and to study the properties of the eigenstates Inverse Participation Ratios (see definition below in section UH which 
have proven to be very appropriate order parameters of Anderson transitions for electronic models (see the review Q) 
The paper is organized as follows. In section UH we introduce the phonon model and the notations for the useful 
observables. Our numerical exact diagonalization results are described in the remaining sections. In section Mil 
concerning the one-dimensional case, the finite-size scaling analysis in the low-frequency region is in agreement with 
the power-law divergence £(w) oc l/uS 2 of the correlation length near zero-frequency [22],l23|]. In section ITVl concerning 
the two-dimensional case, the finite-size scaling analysis in the low-frequency region is in agreement with the essential- 
singularity divergence In £(w) oc 1/u 2 of the correlation length near zero- frequency [23] . Finally in sectionlVlconcerning 
the three-dimensional case, the finite-size scaling analysis around the finite critical frequency is compatible with the 
universality class of the Anderson transition for the Anderson tight-binding electron model in dimension d = 3. Our 
conclusions are summarized in section IVT1 
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II. MODEL AND OBSERVABLES 



A. Scalar phonon problem in a crystal of random masses 

In dimension d = 1, 2, 3, we consider L d random masses mp whose positions f = (m , ..rid) at rest form an hypercubic 
lattice (rii — 1, 2, .., L). These masses are coupled by spring constants Kpp i = 1 if r and r ' are neighbors on the 
hypercubic lattice, so that each mass in the bulk has (2d) neighbors. We consider the following harmonic Hamiltonian 
for the scalar displacements up(t) 

r r,r ' 

The scalar assumption is very standard to simplify the analysis (27j and means physically that longitudinal and 
transverse vibrations are decoupled. Equivalently, the model can be defined by the equations of motion 

mpilp = — Kftf i {up — up /) (2) 

r ' 

To avoid the free-motion of the center-of-mass of the system, we have chosen to consider the fixed boundary conditions 
u = on the lines n.$ = and rij = L + 1 surrounding the hypercube. 

Finally, to avoid the peculiarities introduced by a binary distribution of the disorder [27l . [28[ , we have chosen to 
consider the continuous flat distribution for the random masses mp. 

P(m) = =0 (1 < m < 1 + W) (3) 

so that W represents the disorder strength. The numerical results presented in this paper correspond to the two cases 
W = 1 and W = 20. 



B. Eigenmodes analysis 



Since the equations of motion of Eq. [2] are linear, the dynamics can be analyzed via the eigenmodes of oscillations 
in e lut : the eigenvalues ujp and the associated eigenmodes a p (r) satisfy 



irWpCip(r) = E K P,r ' {a p {r) - a p (f ')) 



(4) 



It is more convenient to perform the similarity transformation }21| 



to reduce the problem to the diagonalization of a symmetric operator 



r ' v 



(5) 



(6) 



As stressed in this form is analog to an Anderson tight-binding model with on-site energies e(r) = , 



and hoppings ■— ^^=, but as a consequence of correlations through the random masses, different physical properties 

can occur. In particular, the eigenvalues are positive E p — ujz > in the phonon problem, whereas E = is the 
center of the band in usual Anderson tight-binding models. It is convenient to work with the orthogonal basis 4> p of 
eigenvectors of Eq. [S] normalized to 



< <t> p \4> P >= E^p( f ) = 1 
f 

This means that the phonons eigenmodes of Eq. [5] are normalized according to 

1 = ^2m F a 2 p (r) 



(7) 



(8) 
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C. Inverse Participation Ratios (I.P.R.) 

To characterize the localization properties of the phonon eigenmodes a p (r) introduced above in Eq. [4[ we consider 
the Inverse Participation Ratio (I.P.R.) 

Y^L^J^W (g) 

that represents an order parameter for Anderson localization transition Q : at a given frequency ui, localized eigen- 
states correspond to a finite value in the limit L —¥ +00 

Y^ oc (uj,L) cx F 2 (u;,oo)>0 (10) 

L— ¥+00 

whereas delocalized states correspond to the following power-law decay 

Y**°°{u,,L) cx -L (11) 

L— > + oo L 

Note that for phonons, the standard definition of the I.P.R. Y% is Eq. [5] in terms of the modes a p (r) that are 
normalized with Eq. [SJ whereas in electronic tight-binding models, the I.P.R. are defined in terms of the orthogonal 
basis <fip Q normalized with Eq. [Jj Since the random masses appearing in the normalization of Eq. [8] remain bounded 
(see the distribution of Eq. [31), we believe that the choice of the a p or of the <\> v to compute the I.P.R. should of 
course affect their precise numerical values, but should not change their scaling properties with the system size L. 
In particular, the localization and derealization criterions of Eqs ITTJ1 and [TT1 should give the same results for the two 
definitions. In the following, all numerical results correspond to the definition of Eq. [9] 



D. Average over the disordered samples of a given size 

In practice, for each size L in dimension d, we generate a certain number ns(L) of disordered samples containing L d 
random masses. The exact diagonalization of each sample a (via the standard NAG diagonalization routine F02FAF 
that computes all the eigenvalues and all the eigenvectors of a real symmetric matrix) yields the L d eigenmodes 

that are ordered by their frequency (dp** in ascending order < uj[ a ^ < lo^ < ... < . For each index p — 1, .., L d , 
we have computed the typical frequency 



ut yp (L) = e lnJ " a) (12) 



and the corresponding typical I.P.R. 



K iw (L)Ee lnF '°( i > (13) 



where A denotes the average of the observable A over the disordered samples (a). The integrated density of states is 
then obtained as 

JVH = i^(^(L)< W ) (14) 
P =i 

with the boundaries values N(0) = and A(+oo) = 1. The parametric plot (oj p vp (L),Y* vp (L)) with p = l,2,..L d 
allows to obtain the behavior of the I.P.R. l^ ti/p (a;, L) as a function of the frequency ui and of the length L. 

III. LOCALIZATION PROPERTIES OF PHONONS IN DIMENSION d = 1 

In this section, we present our numerical results obtained in dimension d = 1. for the following sizes L and the 
corresponding number n s (L) of disordered samples 

L = 10 2 ,2.10 2 ,5.10 2 ,10 3 ,2.10 3 ,3.10 3 ,4.10 3 ,5.10 3 
n„(L) = 2.10 7 ,47.10 5 ,43.10 4 ,4.10 4 ,3.10 3 ,10 3 ,350, 150 (15) 
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A. Density of states 




FIG. 1: (Color on line) Integrated density of states N(u) of Eq. [T4]for phonons in d — 1 (a) N(u) for various sizes 100 < L < 
5000 and two disorder strengths W = 1 and W = 20 (b) same data in log-log scales to display the low-frequency behavior of 
Eq. [16] 

We show on Fig. [T] (a) the integrated density of states N(w) of Eq. [14] for two disorder strengths W = 1 and 
W = 20. As shown in log- log scale on Fig. Q] (b), we find the linear behavior already present in the pure case 

N(uj) ~ C(W)lo (16) 

and the disorder strength W is only present in the numerical prefactor C(W). We have also checked that the lowest 
frequency mode scales as uji(L) oc 1/L. 



B. Typical Inverse Participation Ratio Y 2 vp ((jJ,L) 



To analyse the localization properties of eigenstates, we show on Fig. [2] the typical Inverse Participation Ratio 
Y 2 vp (uj, L) of Eq. [T2] as a function of the frequency u> for various sizes L. In the high-frequency domain where the 
data of all sizes collapse, the eigenstates are localized. In the low-frequency domain where all sizes give different 
results, as shown more clearly in log- log scale on Fig. [5] (b), eigenstates are delocalized on the whole disordered 
sample. We find moreover that the data for the very different disorder strengths W = 1 and W = 20 merge in the 
low-frequency region for each size L : this means that the lowest frequencies eigenstates are delocalized in the same 
way independently of the disorder strength W. 



C. Finite-size scaling analysis of the low-frequency modes 



We show on Fig. [3] that our data are compatible with the following finite-size scaling for the low-energy modes 



V*yp 
1 1 



(uj,L)~jF d=1 (^L 1/u u?j with^ = 2 (17) 
This means that to each frequency lu, one can associate a correlation length diverging as the power-law 

£M « A (18) 

in agreement with 22, 23] . An eigenmode of fixed frequency u> will be delocalized on samples of small lengths L <C £(w) 
with an I.P.R. of order 

Y?>(u>,L) * ^=i(0) (19) 
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FIG. 2: (Color on line) Typical Inverse Participation Ratio Y2(lo,L) as a function of the frequency uj in d — 1 for various sizes 
100 < L < 5000 (a) Y* vp (uj,L) as a function of uj for the disorder strength W = 1 (b) In ^"'(w, L) as a function of lno; for 
two disorder strengths W = 1 and W = 20 




_4 _3 _2 _1 1 2 3 4 5 -5 -4 -3 -2 -1 1 2 3 4 



FIG. 3: (Color on line) Finite-size scaling analysis of the typical I.P.R. Y 2 tap (aj,L) of the low-frequency eigenmodes in d = 1 
according to Eq. 1171 : the rescaled variable y = LY^iuj, L) is plotted as a function of the reduced variable x = L^^uj with 
the value v = 2 (a) data collapse in log-log scale for W = 1 (b) data collapse in log-log scale for = 20 



but will be localized on samples of large lengths L 3> with an I.P.R. of order (using Fd = i(a;) oc at large x) 

Y* w {uj,L) ~ w 2 = 7^t (20) 

The conclusion is thus that for any fixed frequency uj, the corresponding eigenmodes will become localized in the limit 
L — > +oo [22|, [23|. However if one is interested into the set of eigenstates of a sample of a given size L, the conclusion 
is that frequencies u > L -1 / 2 correspond to localized eigenmodes, whereas a certain number Mdeioc{L) of eigenmodes 
corresponding to frequencies uj < L^ 1 / 2 are delocalized. From the linear behavior in w of the integrated density of 
states of Eq. [THl one obtains that the fraction of delocalized states scales as the pseudo-critical value ui*{L) ~ £~ 1//2 

Afdeloc(L) g^),^-!^ (21) 
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So the number of delocalized eigenstates in a sample of size L grows as 

N d eioc(L) oc L 1 / 2 (22) 

L— s-+oo 

in agreement with the interpretation given in [22| . Note that this property of phonons in dimension d = 1 is very 
different from the Anderson electronic problem, where the whole set of eigenstates of a given sample become localized 
at large sizes. Physically, this difference is essential if one considers the dynamical properties, since the dynamics 
can be expanded on the basis of eigenmodes : in the Anderson electronic problem, the localization of the whole set 
of eigenfunctions imply the exponential localization for the dynamical problem starting from any localized initial 
condition, whereas in the phonon case, the presence of these low-frequency delocalized modes for any size prevents 
the exponential localization in the dynamics. 



IV. LOCALIZATION PROPERTIES OF PHONONS IN DIMENSION d = 2 



In this section, we present our numerical results obtained in dimension d = 2. for the following sizes L and the 
corresponding number n s (L) of disordered samples 

L = 10,20,30,40,50,60,70,80 
n s (L) = 2.10 7 ,7.10 5 ,77.10 3 , 7500, 1800,500,250,200 (23) 



A. Density of states 




FIG. 4: (Color on line) Integrated density of states N(tj) for phonons in d — 2 (a) N(u) for various sizes 10 < L < 80 and two 
disorder strengths W = 1 and W = 20 (b) same data in log-log scales to display the low-frequency behavior of Eq. 1241 

We show on Fig. [4] (a) the integrated density of states N(uj) of Eq. [14] for two disorder strengths W = 1 and 
W = 20. As shown in log- log scale on Fig. |4] (b), we find the same behavior as in the pure case 

N(u) cx u? (24) 

and the disorder strength W is only present in the numerical prefactor. We have also checked that the lowest frequency 
mode scales as u>i(L) ocl/L, 



B. Typical Inverse Participation Ratio Y 2 vp (u),L) 

We show on Fig. [5] the typical Inverse Participation Ratio F 2 tl,p (w,L) of Eq. Q2]as a function of the frequency u> 
for various sizes L. In the high-frequency domain where the data of all sizes collapse, the eigenstates are localized. In 
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FIG. 5: (Color on line) Typical IPR Y 2 vp (aj, L) as a function of the frequency co in d = 2 for all eigenmodes of samples of sizes 
20 < L < 80 (a) Y* yp (uj,L) as a function of uj for W = 1 (b) In F 2 t!/P (aj, L) as a function of In a; for two disorder strengths 
W = 1 and W = 20 



the low- frequency domain where all sizes give different results, as shown more clearly in log- log scale on Fig. [3] (b), 
eigenstates are delocalized on the whole disordered sample. As in dimension d = 1, we find moreover that the data 
for the two disorder strengths W = 1 and W = 20 merge in the low-frequency region for each size L : this means that 
the lowest frequencies eigenstates are delocalized in the same way independently of the disorder strength. 



C. Finite-size scaling analysis of the low-frequency modes 




In ( 0) (In L) 



1/(1 



In ( a (In L) 



FIG. 6: (Color on line) Finite-size scaling analysis of the low-frequency modes in d — 2 for 30 < L < 80 according to Eq. [25] : 
y = L 2 Y2 VP (lo, L) as a function of z = (lnL) 1//M cj with the value (i = 2 (a) data collapse in log-log coordinates for W = 1 (b) 
data collapse in log-log coordinates for W — 20 

We show on Fig. [5] that our data are compatible with the following finite-size scaling for the low-energy modes 

Y* vp {u,L) ~ ^F d=2 ((lnL) 1 /^) with ^ = 2 (25) 



8 



This means that to each frequency ui, one can associate a correlation length diverging as the following essential 
singularity 

In CM oc i (26) 

in agreement with (23j . As in dimension d — 1, the conclusion is thus that any fixed frequency mode u> becomes 
localized in the limit L — > +oo [23$ . 

However if one is interested into the set of eigenstates of a sample of a given size L 2 , the conclusion is that frequencies 
uj > (InL) -1 / 2 are localized, whereas a certain number Afdeioc(L) of eigenstates with frequencies w < (InL) -1 / 2 are 
delocalized. From the behavior in u of the integrated density of states of Eq. [24l one obtains that the fraction of 
delocalized states scales as the square of the pseudo-critical value ~ (InL) -1 / 2 

M del oc{L) oc( ^ (£)) 2^ (lnL) -l (27) 

So the number of delocalized eigenstates in a sample of size L 2 grows as 

L 2 

Ndeloc(L) cx — (28) 

l^+oo In L 



V. LOCALIZATION-DELOCALIZATION TRANSITION OF PHONONS IN DIMENSION d = 3 



In this section, we present our numerical results obtained in dimension d — 3 for the following sizes L and the 
corresponding number n s {L) of disordered samples 

L =8,10,12,14,16,18 
n s {L) = 35.10 4 ,45.10 3 , 5500, 1500, 300, 300 (29) 



A. Typical Inverse Participation Ratio Y 2 vp (lu, L) 




FIG. 7: (Color on line) Typical IPR F 2 ' ap (cj, L) as a function of the typical frequency in d = 3 for all eigenmodes of samples of 
sizes 8 < L < 18 for W = 20 (a) In Y^iui, L) as a function of lno; : the collapse in the high frequency region corresponds to 
localized states (b) ln(L 3 Y 2 tap (u;, L)) as a function of In a; : the collapse in the low frequency region corresponds to delocalized 
states 

We show our data for the typical I.P.R. y 2 * yp (tj, L) on Fig. [7] (a) : in the high-frequency part of the spectrum, the 
data collapse for the various sizes L corresponds to localized states with a finite value Y 2 vp (lj,oo) > (see Eq. ITU]) . 
We show on Fig. [7] (b) the same data after the appropriate rescaling L 3 F 2 * yp (ti;, L) to detect the delocalized states 
(see Eq. ITTj) : the data collapse in the low-frequency part of the spectrum corresponds to delocalized states. 
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B. Finite-size scaling analysis of the localization transition 




FIG. 8: (Color on line) Finite-size scaling analysis of the localization transition in d — 3 for the disorder strength W = 20 (a) 
\n(L D ( T) Y2(uJ, L)) as a function of In a; for the sizes 8 < L < 18 with the value D{2) = 1.3 : the crossing determines the critical 
point ln.Wc(W = 20) ~ —0.1 (b) ln(L D ^Y2(u>, L)) as a function of ln(o; — ui c ) + (l/^)lnL for v = 1.57 : the data collapse is 
satisfactory. 

In dimension d = 3, one expects that there exists a localization-delocalization transition at some finite frequency 
ujc > [23]. The I.P.R. is then expected to follow the following finite-size scaling 



V, L) ~ j^F d=3 (LV"( W - Wc )) (30) 

The exponent D(2) governs the power-law decay of the I.P.R. exactly at criticality 

^(o^Ljoc-^ (31) 

For the transition of the Anderson tight-binding electronic model in d = 3, it has been measured numerically (see @ 
and references therein) 

DAnderson{2) — 1.3 (32) 

As shown on Fig. [8] (a), if we rescale our data using this value, we obtain that the curves L d ^Y2 VP (uj, L) for various 
L cross around the value lnw c (VF = 20) ~ — 0.1 corresponding to 

uj c {W = 20) ~ 0.9 (33) 

The integrated density of states at this value is around N(uj c (W — 20)) ~ 0.66 (data not shown), so that the critical 
point is sufficiently inside the spectrum to have enough localized states and delocalized states on both sides (this is 
not the case for any value of the disorder strength as explained below in section IV C|) . 

In addition, if we now rescale our data in terms of the reduced variable L 1 / I/ (a; — uj c ) with the value of the correlation 
exponent that has been measured numerically for the Anderson tight-binding electronic model in d = 3 (see 0] and 
references therein) 

V Anderson — 1.57 (34) 



we obtain a good data collapse as shown on Fig. [5] (b). Our conclusion is thus that the localization transition of 
phonons in d = 3 is governed by the same universality class as the Anderson tight-binding electronic model in d = 3. 



10 




In CO In co 



FIG. 9: (Color on line) Typical IPR as a function of the frequency to in d = 3 for all eigenmodes of samples of sizes 8 < L < 18 
for the disorder strength W — 1 (a) In Y 2 yp as a function of In lj : even in the region of the highest frequencies, there is no 
collapse that would correspond to localized states (see the difference with the data corresponding to W — 20 on Fig. [7] (a)) 
(b) ln(L 3 Y^ HP ) as a function of lnw : the collapse on the whole frequency region means that all states are delocalized (see the 
difference with the data corresponding to W — 20 on Fig. [3(b)) 

C. On the importance to consider strong enough disorder W to observe the transition 

Up to now, we have described our results in dimension d = 3 for the disorder strength W — 20, for which we have 
found a clear localization transition at some finite frequency uj c (W = 20) (Eq. I33|) well inside the spectrum. However, 
as in dimensions d = 1 and d = 2, we have also studied the disorder strength W = 1 : the corresponding data for the 
typical I.P.R. Y 2 typ (uj, L) shown on Fig. [9] indicate that here all states of the spectrum are actually delocalized. This 
does not mean that there is no critical frequency ui c (W = 1), but that this critical value is not accessible, because it 
is higher that the maximal frequency uj max of the spectrum where the density of states reaches its asymptotic value 
N(u} max ) = 1. Our conclusion is thus that to observe the localization transition, one should consider sufficiently 
strong disorder W to ensure that the corresponding critical value belongs to the spectrum 

u c (W) <uj max (W) (35) 

We have found that this condition is satisfied for W = 20, but is not satisfied for W = 1. Changing the value of W 
allows to move the critical value ui c (W) inside the spectrum. 

VI. CONCLUSION 

To characterize the localization properties of cigcnstates for phonons in the presence of random masses in dimension 
d = 1, 2, 3, we have studied numerically the behavior of the typical Inverse Participation Ratio l^O^; L) as a function 
of the frequency ui and of the linear length L of the disordered samples. 

In dimensions d — 1 and d = 2, we have found that the low- frequency part uj — > of the spectrum satisfies 
the following finite-size scaling LY 2 (u},L) = Fd=\(L x ^ 2 iS) in dimension d — 1 and L 2 Y 2 (u),L) = Fd= 2 ((\n L) 1 / 2 ^) in 
dimension d = 2. We have moreover explained that the loose statement "all eigenstates are localized in dimensions 
d = 1,2" should be stated with some care for phonons : it is true that an eigenstate of any fixed frequency w becomes 
localized in the limit L — > +oo, but one should also be aware that a given disordered sample of fixed length L contains 
a certain number Nd e i oc (L) of delocalized states growing as Ndeioc(L) ~ L 1 / 2 in d = 1 and as Ndeioc(L) ~ L 2 / (InL) in 
d = 2. These low-frequency delocalized modes are expected to play a major role in the dynamical properties on large 
distances, and in particular in the heat transport problem where the disordered sample is connected to heat baths at 
the boundaries (see for instance (27j). 

In dimension d = 3, for strong enough disorder strength W (W — 20 in our case), we have found a very clear 
localization-delocalization transition at some finite critical frequency w c (W) > 0. We have shown that our data are 
compatible with the finite-size scaling L d ^Y 2 {uj,L) = F d=3 (i 1/l/ ( 

lu — uj c )) with the values D(2) ~ 1.3 and v ~ 1.57 
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corresponding to the universality class of the localization transition for the Anderson tight-binding electronic model 
in dimension d = 3. We have also found that for too small disorder strength (namely W = 1 in our case) the 
critical point uj c (W) can be higher that the maximal frequency of the spectrum, so that all eigenstates are actually 
delocalized. The choice of too small disorder strengths seems to be the reason why localized states were found only 
very near band-edges in previous numerical studies [2J-|27| . 
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